Cooling of a suspended nanowire by an AC Josephson current flow 
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Abstract 

We consider a nanoelectromechanical Josephson junction, where a suspended nanowire serves as a superconducting weak 
link, and show that an applied DC bias voltage can result in suppression of the flexural vibrations of the wire. This cooling 
effect is achieved through the transfer of vibronic energy quanta first to voltage driven Andreev states and then to extended 
quasiparticle electronic states. Our analysis, which is performed for a nanowire in the form of a metallic carbon nanotube and 
in the framework of the density matrix formalism, shows that such self-cooling is possible down to the ground state of the 
flexural vibration mode of the nanowire. 

PACS numbers: 73.23.-b, 74.45.+C, 85.85.+j 



i 

C 

o 
o 



(N 
> 

o 

(N 

(N 
O 
O 



X 



Nanoelectromechanical systems (NEMS) have over the 
last two decades been a very active field of both funda- 
mental and applied research. With typical dimensions on 
the nanometer scale NEMS combine electronic and me- 
chanical degrees of freedom for novel applications. These 
include small mechanical resonators for ultrasensitive 
mass spectroscopy {1:!, 2] and position sensing with close 
to quantum limited resolution 0, 0] ■ The prospect of a 
controlled fabrication of NEMS devices whose physics is 
ultimately governed by the laws of quantum mechanics 
is a major driving force in the field today. 

Typically, nanoelectromechanical systems comprise a 
mechanical resonator coupled to an electronic system 
used for both actuation and detection. Due to the high 
mechanical vibration frequencies and the exceptionally 
high quality factors recently achieved, these systems al- 
low for very low energy dissipation and extreme sensi- 
tivity to external stimuli. To make full use of their po- 
tential, much research has recently focused on the pos- 
sibility to cool the mechanical subsystem in NEMS to 
its vibrational ground state. Several theoretical papers 
have invoked side-band cooling or dynamical back-action 
in order to predict electromechanically induced cooling 
of mechanical resonators to the ground state 5-8] . The 
best experimental result so far, corresponding to a vibron 
occupation number of 3.8, is due to Rocheleau et al. 
They used an external electromagnetic field as an energy 
transducer for side-band cooling of the mechanical res- 
onator. Similar results have been achieved by cooling 
with dynamical back-action [l(| ■ 

Here, we propose a new cooling method, which utilizes 
the unique electronic properties of a Josephson junction. 
We consider an SNS junction which is composed of a 
nanowire in the form of a metallic carbon nanotube sus- 
pended between two superconducting leads as in Fig. [T] 
Below we demonstrate that such a nanomechanical weak 
link has the capacity to self-cool, if biased by a DC volt- 




FIG. 1: (Color online) Schematic diagram of the Josephson 
junction discussed in the text. A suspended carbon nanotube 
of length L forms a weak link between two superconducting 
leads. An external magnetic field H gives rise to a Lorentz 
force that couples flexural vibrations of the nanowire to the 
electrical current that is carried by Andreev states through 
the junction if it is biased by a voltage V. Cooling is achieved 
by the absorption of vibration-energy quanta as the electronic 
subsystem is excited across a gap in the Andreev spectrum 
(see Fig.©. 



age V, through the transfer of energy from the flexural 
vibrations of the nanowire to a bath of electronic quasi- 
particle excitations. Recently, it has been experimentally 
shown that for systems similar to the one considered here 
it is possible to cool nanomechanical objects below the 
ambient temperature through the transfer of mechanical 
energy to electronic quasiparticle states [ill E3 • Similar 
to these works, the cooling process to be proposed is an 
inherent property of the system, which makes the sug- 
gested mechanism a promising candidate for truly quan- 
tum mechanical manipulation of mechanical resonators. 

The electromechanical coupling required for our cool- 
ing mechanism can be achieved by applying an exter- 
nal magnetic field perpendicular to the length of the 
nanowire. Displacement of the nanowire then couples 
the mechanical and electronic degrees of freedom through 
the Lorentz force. For an SNS junction that is short com- 
pared to the coherence length in the superconductors, the 
AC Josephson dynamics induced by the weak DC bias 
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voltage can be expressed in terms of a pair of Andreev 
states, which is periodically created and destroyed in the 
junction. Between creation and destruction the Andreev 
levels evolve adiabatically within the energy gap of the 
superconductors; the level spacing first shrinks, reaches 
a minimum and is then restored to its initial value as the 
Andreev states are dissolved in the continuum quasipar- 
ticle spectrum (see Fig. [2]). Being thermally populated 
when created, the population of the Andreev levels be- 
come "over-cooled" as they evolve towards the middle of 
the gap. Since the Andreev states are coupled to the 
nanowire vibrations, they can absorb a quantum of vibra- 
tion and approach an equilibrium distribution through 
an electronic excitation from the lower to the upper An- 
dreev level. The absorbed energy is transferred to the 
electronic quasiparticle continuum (and therefore trans- 
ported away from the junction) when the Andreev states 
are dissolved. Cooling of the nanowire proceeds by re- 
peated absorptions of vibrational quanta in subsequent 
cycles of the temporal evolution of the Andreev levels. 

To discuss the cooling mechanism quantitatively we 
consider the model Hamiltonian, 
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H(t) = £{4>{t))a z + A VR8in((t>(t)/2)a x + 
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Here, the first two terms describe the adiabatic dynam- 
ics (h(f> <C Ao) of the two level system formed by the 
Andreev states that carry the Josephson supercurrent 
through the (normal) nanowire HI, 14 1 ; a a are Pauli ma- 
trices, Ao is the order parameter in the superconducting 
leads, R is the normal state electronic reflection proba- 
bility, and £(<fi(t)) = Ao cos(0(t)/2) is the energy of the 
Andreev states for the completely transparent junction 
with <p(t) the phase difference between the two supercon- 
ductors. The latter depends on the bias voltage according 
to the Josephson relation cb — 2eV/ H. The third term in 
the Hamiltonian models the vibrating nanowire as a har- 
monic oscillator, taking only the fundamental bending 
mode into account; the operator w [b] creates [annihi- 
lates] one quantum huj of vibrational energy. The last 
term in (JT|) describes the electromechanical coupling of 
the current-carrying Andreev states to the motion of the 
nanowire through the Lorentz force. Here, e is the elec- 
tronic charge, L is the length of the wire, yo is the zero- 
point oscillation amplitude, H is the magnetic field and 
2e/h(d£ I d<j))d z is the current operator in the nanowire 

In the adiabatic regime, i.e. when eV < eV c — 4i?Ao, 
it is convenient to use a basis set formed by the states 
ip±((j>(t)) with corresponding energies E±(t) = ±E(t), 
where E(t) = A (l - Dsin(0(i)/2)) 1 / 2 and D = 1 - R 
[l6j |. In this basis, the Hamiltonian (p} reads, 

H e ff(t) = E{t)f z + tujtfb + A $(tf + b)f x . (2) 
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FIG. 2: a. Temporal evolution of voltage-driven Andreev 
levels (full lines) in the Josephson junction of Fig. [T] Here, 
j±, n) are states of the junction; |±) denotes the upper/lower 
Andreev level and |n) is the quantum state of the nanowire. 
Transitions between | — , n) and | + ,n — 1), through the emis- 
sion/absorption of one vibrational quantum, can be modeled 
by a scattering matrix S. The transition probabilities depend 
on the state of the nanowire through |^2(w)| 2 = 1 — |k2(w)| 2 . 
After one period, TV, the partially filled Andreev levels join 
the continuum, creating quasiparticle excitations (represented 
by dashed arrows) and the electronic states are reset (filled 
and empty circles). Cooling of the nanowire occurs because 
absorption of vibrational quanta is more probable than emis- 
sion at the ambient temperature, b. Corresponding energy 
flow diagram. In each cycle the voltage-driven Andreev levels 
(AL) move some energy E from the mechanical resonator to 
a quasiparticle bath in thermal contact with a heat bath at 
temperature T. In this process the battery supplies an energy 
W . Only if the electronic and mechanical sub-systems are 
coupled can E and W be non-zero. The resonator is coupled 
by dissipation to the same heat bath from which it receives 
an energy E^{= E in the steady state). "Cooling" ["heating"] 
refers to decreasing [increasing] the energy separation between 
the Andreev levels. 



Here, the Pauli matrices Ti span the space formed by the 
states ii±(4>(t)) (see Fig. [2]). In the last term of ((2]), $ = 
2LHiryo / $o is the dimensionless magnetic flux threading 
the area swept by the vibrating nanowire; $0 — h/2e is 
the magnetic flux quantum. 

One notes that when cf> = n, the Andreev states 
'4 > ±{4 > = n )> with energies E±(cf> = n) = ±AqR 1 / 2 , are 
symmetric/antisymmetric superpositions of states carry- 
ing current in opposite directions. Transitions between 
the Andreev states — induced by the nanowire as it 
vibrates in a transverse magnetic field — are therefore 
by far most probable when ip(t = to) = n. Below 
we will therefore consider the resonant situation when 
Huj = 2A0R 1 ' 2 [17|, which is the optimum condition for 
the proposed cooling mechanism. To proceed further we 
apply the rotating wave approximation, after which the 
Hamiltonian reads 



E(t) - hu/2 
A <S>ti 



A $6 

-E{t) + huj/2 



(3) 



Within the conditions for adiabaticity outlined above, 
we can now analyze the probability for mechanically in- 
duced transitions between the Andreev levels as a func- 
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tion of the coupling strength. As such transitions occur 
only in the vicinity of to, one can use the parabolic ex- 
pansion of E(t) around this point, which leads to the 
following time-dependent Schrodinger equation, 



id T c n (T) = (t 2 ct 2 + Ty/na x ) c„(r) 
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Here, c„(r) = (c +>n _i(r), c^^(t)) t where c±,„(t) is the 
probability amplitude for finding the system in the up- 
per/lower electronic branch with the harmonic oscillator 
in the state n. 

Within the adiabatic regime, where V < V c , the 
coupling between the Andreev levels is relatively weak; 
T = Ao$/(fiu;) < 1 (see below). This allows us to 
analyze (0} by treating T as a perturbation to the so- 
lution |c_ !rl (r)| 2 = 1, where the system starts in the 
lower electronic branch at time — 5t = —St^/h) 1 / 3 and 
remains there until +St. Solving (jj) for the prob- 
ability |c +! „_i((5r)| 2 to find the system in the upper 
branch after passing through the resonance, we find that 
|c+. n _i(<5r)| 2 ~ imT 2 which is consistent with our nu- 
merical analysis. 

In order to model the full evolution of the coupled elec- 
tromechanical system we evaluate the total density ma- 
trix p of the system over one period. Under the assump- 
tion that the external temperature, T, is much smaller 
than the superconducting gap, the system will initially 
be found in the lower electronic branch with n quanta of 
the nanowire vibrations excited with probability P™, 
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During the adiabatic evolution (no coupling between 
branches) the system interacts with the external heat 
bath and the rate of change of the density matrix is given 
by the expression, 



dp(t)/dt-- 



-i[je eff (t),p{t)]/* + 'yHW))/2 

-(1 + n B ) (i> r bp + pPb - 2bptf 



(5) 



ffip + pbP - 2ppPj , (6) 



where the collision integral £(/5) models the interaction of 
the mechanical subsystem with the environment. Here, 
n B = (exp(/3t>w) - I)" 1 with /3 = (ksT)- 1 while 7 = 
lo/Q is the damping rate of the nanowire vibrations with 
Q the quality factor. Considering the coupling to the 
environment to be small (Q ~10 5 [lj|) we can proceed by 
treating the interaction of the system with the external 
environment as a perturbation. 



The large separation of the energy scales ensures that 
transitions between the Andreev levels is almost always 
negligible. It is only during a short time of order St ~ 
(Sw$/A ) 1/2 /i/(eVO < nh/(eV) around t that the two 
branches interact, which permits us to describe the evolu- 
tion of p through the resonance, p(t +5t) = Sp(t — St)^^, 
by the unitary scattering matrix S, 



S = 



Ki(n) 

Vn+l 



v\ (ft) 



Vn+l 



(7) 



In (|7|), h — b^b is the vibron number operator and the 
subscripts 1, 2 refer to the top/bottom Andreev level re- 
spectively. The coefficient Ki(n) [vi(n)] is the probability 
amplitude for the system to stay in [scatter out of] the 
initial state i, which depends on the state of the oscilla- 
tor as outlined above {\ni{n)\ 2 + \i>i(n)\ 2 = 1). As such, 
I v 2{n)\ 2 is the probability of the system, initially in the 
lower electronic branch with n vibrons excited, to scat- 
ter into the upper electronic branch through the absorp- 
tion of a vibron. It thus corresponds to |c+ !rt _i(^r)| 2 in 
(dJ. We conclude that |^2 | 2 — TrnT 2 and note that 
|^2(n)| 2 = |^i (n — 1)| 2 , from the symmetry of (j4|). 

After one period the bound Andreev states merge with 
the continuum of itinerant states. At this time the quasi- 
particle carried by the Andreev states is released into the 
continuum and the initial conditions of the Andreev level 
populations is reset (see Fig. [2]) 20]. The mechanical 
system will, however, not return to the initial vibrational 
state. Instead we must find its final state from the den- 
sity matrix for the mechanical system after one period, 
Pmech' obtained by tracing out the electronic degrees of 
freedom in the total density matrix, 

/SLcfc = IMP 7 ) = Tr ei (Sp m & + 



i^-sz(p m )& + l^-Hsrs^)) . in) 



Equation © describes the evolution of the mechani- 
cal density matrix over one period. Without coupling to 
the environment, 7 = 0, this corresponds to a decay of 
the mechanical subsystem, since for each period there is 
a probability |i/ 2 («)| 2 that the state |— ,n) scatters into 
the state |+, n — 1). As the opposite process is forbidden 
if initially only the lower electronic branch is populated, 
the mechanical subsystem would thus approach the vi- 
brational ground state at a rate, which depends on the 
strength of the coupling. Including effects of the environ- 
ment, however, the mechanical subsystem will not reach 
its ground state as the external damping drives the sys- 
tem towards thermal equilibrium. In Fig. [3] we plot the 
evolution — over many periods starting from its initial 
thermal distribution — of the probability P„ for n quanta 
of mechanical nanowire vibrations to be excited using re- 
alistic experimental parameters. As can clearly be seen, 
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should act as a self-cooling device as the "over-cooled" 
Andreev states can readily absorb energy from the me- 
chanical subsystem given that sufficient coupling between 
the two can be achieved. The corresponding energy up- 
take of the electronic subsystem is later released into the 
quasiparticle continuum, leading to an effective cooling 
of the nanomechanical resonator. 
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FIG. 3: (Color online) Distributions of the probabilities P„ 
for n quanta of the lowest flexural vibration mode of the sus- 
pended nanowire in Fig. [T]to be excited after N of the cool- 
ing cycles described in Fig. [2] Initially P n is thermally dis- 
tributed, P n = exp(— /3hum)(l — exp(— /3ftcj)), with fcsT = 
5tku. The plots were calculated for hu) = 10 -6 eV, Ao = lOftcj, 
yo =20 pm, L =100 nm, H =1T and T v ~20ns. The inset 
shows | V2 (n) | 2 for the same parameters. 



the scattering process acts to lower the number of excited 
quanta, effectively cooling the nanowire down to a final 
average vibron population (n) ~ 0.1. This is confirmed 
by a perturbative calculation with T <C 1, which shows 
that to first order in the damping rate, 7, the stationary 
population of the first excited state is 



1 hu 



1/3 



1 



ophu 1 



0.1, (9) 



which is in accordance with the distribution in Fig. [3] 
Here we note that including a small initial population 
of the upper Andreev branch, A" < 1, only changes the 
stationary vibron population to ~ (n) +0(c™). Thus, in- 
cluding, e.g., a thermal population of the upper Andreev 
branch only increases the stationary vibron population by 
a factor c™ ~ exp(— 2(3A ) ~ 0.02, and does not change 
the qualitative results of this paper. 

To conclude, we have shown that quantum mechani- 
cal cooling of a nanomechanical resonator acting as the 
weak link in a Josephson junction is possible. In par- 
ticular, we have considered the example of a suspended 
carbon nanotube, where the unique combination of high 
resonance frequencies and high mechanical quality factors 
together with a high transparency for electrons combine 
to give sufficiently strong coupling for efficient cooling. 
Using realistic experimental parameters we have shown 
that for a short suspended nanotube, a stationary distri- 
bution of mechanical vibrational quanta corresponding 
to an average occupation number as low as (n) = 0.1 
can be achieved. This is truly in the quantum mechani- 
cal regime. Furthermore, the suggested mechanism does 
not rely on any external electromagnetic fields to stim- 
ulate the cooling process. Rather, the proposed system 



* Electronic address: gustav.sonne@physics.gu.se 
[1] B. Lassagne et al, Nano Lett. 8, 3735 (2008). 
[2] K. Jensen, K. Kim, and A. Zettl, Nat. Nanotechnol. 3, 
533 (2008). 

[3] M. D. LaHaye et al, Science 304, 74 (2004). 

[4] S. Etaki et al, Nat. Phys. 4, 785 (2008). 

[5] I. Martin et al, Phys. Rev. B 69, 125339 (2004). 

[6] S. Zippilli, G. Morigi, and A. Bachtold, Phys. Rev. Lett. 

102, 096804 (2009). 
[7] I. Wilson-Rae, P. Zoller, and A. Imamoglu, Phys. Rev. 

Lett. 92, 075507 (2004). 
[8] S. H. Ouyang, J. Q. You, and F. Nori, Phys. Rev. B 79, 

075304 (2009). 
[9] T. Rocheleau et al, Nature 463, 72 (2010). 
[10] A. Naik et al, Nature 443, 193 (2006). 
[11] P. J. Koppinen and I. J. Maasilta, Phys. Rev. Lett. 102, 
165502 (2009). 

[12] J. T. Muhonen et al, Appl. Phys. Lett. 94, 073101 
(2009). 

[13] P. F. Bagwell, Phys. Rev. B 46, 12573 (1992). 

[14] C. W. J. Beenakker, Phys. Rev. Lett. 67, 3836 (1991). 

[15] L. Y. Gorelik et al, Phys. Rev. Lett. 75, 1162 (1995). 

[16] In this regime, the volt age- induced Landau-Zener tran- 
sition at <f> = 7r between the adiabatic Andreev levels is 
suppressed, ensuring that without electromechanical cou- 
pling the level population remains constant in time. 

[17] This condition can be achieved by controlling the mag- 
nitude of the order parameter Ao through the ma gnet ic 
field for junctions with small reflection coefficients [la ]. 

[18] J. Kong et al, Phys. Rev. Lett. 87, 106801 (2001). 

[19] A. K. Hiittel et a^Nano Lett. 9, 2547 (2009). 

[20] As explained in ;2l|, the non-equilibrium population of 
the Andreev states that approach the edges of the quasi- 
particle energy spectrum at E = ±Ao (Fig. [2]) cannot, 
due to selection rules, affect the population of the new 
Andreev states, which are created with populations "re- 
set" to the thermal equilibrium and move away from the 
edges. This system will not be heated since (i) the re- 
leased energy is carried away by quasiparticles [22J , which 
escape from the vicinity of the weak link and spread over 
the reservoirs with greatly reduced density and (ii) the 
"escape time", T csc ~ (eV/ Ao) 1//3 Ty, is short compared 
to the period Tv- 

[21] L. Y. Gorelik et al, Phys. Rev. Lett. 81, 2538 (1998). 

[22] E. N. Bratus et al, Phys. Rev. B 55, 12666 (1997). 



4 



